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1. Introduction 

One of the most important lessons that one learns in quantum mechanics courses is to 
choose a proper picture for describing the dynamics of a system. An appropriate picture 
is the one in which the physical properties of the system can be easily revealed and the 
mathematical calculations involved are relatively simple. In general, there are three 
pictures for this aim: the Schrodinger, the Heisenberg, and the Interaction picture pp. 
The Schrodinger picture (SP) is more suitable for studying closed and conservative 
systems. In this picture the only time-dependent quantity of the system is its quantum 
state. Correspondingly, observables of the system are constant in time. However, the 
number of situations in which a quantum system can be considered as a closed one is 
limited and, furthermore, most of the interesting phenomenon in the quantum world 
such as decoherence and optical cooling, to just name two, happen when a quantum 
system couples to its surrounding environment, to an external force or to both of them. 
Hence, the Hamiltonian describing the system is more complicated due to the presence 
of the environment and/or the external force. In many of such situations it is impossible 
to find an exact solution to the Schrodinger equation. However, sometimes, adopting the 
Heisenberg picture (HP) makes it possible to find the time-evolution of the expectation 
values of the observables of the system. In HP the state does not vary with time while 
the observables are time-dependent and their evolution is described by the Heisenberg 
equation of motion. Finally, the last scheme, namely the Interaction picture (IP) one, 
is the most suitable picture when the Hamiltonian of a system can be written as a 
sum of two parts: a time-independent term of which the eigenstates and eigenenergies 
are known, and a usually time-dependent term which influences the dynamics of the 
system. Employing this picture enables one to set the dynamics arising from the time- 
independent Hamiltonian aside and focus on the influence of its time-dependent part on 
the evolution of the system. That is why this picture is frequently employed in quantum 
optics where matter interacts with the radiation field. 

All these three pictures are described in terms of unitary transformations in Hilbert 
space. The aim of this work is to map such transformations in phase space and illustrate 
the simplicity of the IP time-evolution of a quantum system in comparison to its SP 
counterpart. Although SP and HP have been mapped into classical active and passive 
transformations [lj, to the best of our knowledge, there is no such direct comparison for 
the interaction picture in the physics literature. Here we fill this gap and show that by 
taking active local and global transformations one can transform the IP to the SP or 
HP. In this way, it is easier to understand the concept of performing different unitary 
transformations in Hilbert space and its consequences on the evolution of the state of 
the system of interest. 

To this end, we take a single harmonic oscillator as our quantum system, initially 
prepared in an ideal squeezed state [2] and first study its free time-evolution in SP 
and HP. Then by employing the Wigner-Weyl description we calculate the time- 
evolution of its corresponding Wigner function. The basic advantage of the Wigner-Weyl 
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representation [HI SI El E] is that the operators x and p turn into c-numbers making the 
correspondence with canonical transformations possible. The Wigner function of the 
coherent and squeezed states of a quantized harmonic oscillator is a Gaussian function of 
these c— number variables and its corresponding phase space distribution can be pictured 
by circles and ellipses, respectively [7JE1I9]. This is the subject of sections [2] and El In 
section 0] we introduce two types of a time-dependent potential to the system: a linear 
and a quadratic term (in terms of the bosonic creation and annihilation operators). Then 
we take the IP to study the dynamics of the initially squeezed state under the driving 
potentials. We explicitly show that the IP- time-evolution of the state in Hilbert space 
corresponds to a transformation to a local frame that is co-rotating with the Wigner 
function in phase space. Furthermore, we show switching from one picture to another is 
equivalent to performing different active and/or passive transformations in phase space. 



2. Basics of the quantized harmonic oscillator 

A single mode harmonic oscillator with unit mass and frequency uq is described by the 
sum of its kinetic and potential energy 

H = ^(f + x>). (1) 

Here x and p represent quadratures of the quantized harmonic oscillator and are defined 
in terms of the annihilation operator a and the creation operator at as 

x = — -={a + fit), (2a) 
v2 

p = i7! (a ~ at) ' {2b) 

where h~ is set to one. Alternatively, this system may be described in terms of a and a) 
as 

H = w (tfa + , (3) 

where the commutator of the creation and the annihilation operator is [a, at] — 1. Ho 
has the energy eigenvalues E n defined as 

E n = u] > with n= 0,1,2,.-.. (4) 

The energy eigenstates (number states) \n) such that: 

H \n) = E n \n). (5) 

Given H Q , the time-evolution of quadrature operators x and p is given by the Heisenberg 
equation of motion as 

x{t) = x(0) cos(u t) + p(0) sin(o;ot), (6a) 
pit) = — x(0) sm(oj t) + p(0) cos(a;ot). (66) 
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Any superposition of the number states is a solution to the Schrodinger equation 
for a quantized harmonic oscillator described by the Hamiltonian H . There are two 
classes of such superpositions that are of great interest: coherent and squeezed states. 
A coherent state |a) is described as the displaced vacuum state |0) and is generated by 
a canonical transformation of the vacuum state |0) in phase space 

\a)=D(a)\0), (7) 
where D(pt) is the unitary displacement operator defined as [T0~1 [12] 

D(a) = e ( Qat - a * a ). (8) 
The displacement parameter a, is a complex number and may be decomposed to 

a=A=(a + ib), (9) 

Squeezed states are of interest because they provide reduced fluctuations in one 
quadrature, when compared with the coherent state. A squeezed state is described by 
the unitary squeezing operator [21 [131 EES El IE] S(() in a similar fashion. However, 
the difference lies in the fact that S(Q is not linear in a and d* and is defined as 

S(C) = exp [(Cd 2 - C« t2 ) /2] , C = se w , (10) 

where s = |£| is the strength of the squeezing and 9 = arg(£) determines the direction 
along which the squeezing is performed. In phase space the squeezed states correspond 
to Gaussian distributions with unequal widths as opposed to the symmetric distribution 
of the coherent state. [8] 

Throughout this paper, for the sake of simplicity and without the loss of generality, 
we take both the displacement parameter a and the squeezing parameter ( as real 
numbers. Moreover, we assume the squeezed state to be an ideal squeezed state [2]. 

Since we want to find a classical correspondence for the unitary transformations in 
the Interaction picture, we need to map the time-evolution of the ideal squeezed state 
mentioned above into classical phase space and study its time-evolution. This is best 
described with the help of the Wigner function and Weyl transform [H HI [5]. Knowing 
p{t) the state of the quantum system at any time, we can build up the symmetric 
characteristic function [11] Tr[p(t)D(x)] and from there the corresponding SP- Wigner 
function W(x,p;t) as 

W(x,p;t) =ljdtj d V e-^+^ Tr[p(t)D( X )], (11) 

where D(x) is the displacement operator with the parameter x defined as 

X = 4=(-^7 + iO- (12) 
Both integrations are from — oo to oo. 
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3. Schrodinger and Heisenberg pictures vs. passive and active 
transformations 



Before embarking on the classical correspondence of the IP, let us first explicitly show 
for a single-mode harmonic oscillator how SP and HP map to active and passive 
transformations in phase space, respectively. For a harmonic oscillator with the 
Hamiltonian given by ([3]), the unitary time-evolution operator U , which governs the 
evolution of the state, reads 



Uo(t) 



,-lffo* 



(13) 



If the system is initially in an ideal squeezed state, 

p(0) = t){^=)S (s) |0> (0| S (-s) D{-^j=l (14) 

with /i x representing the displacement along the x-axis in phase space and s denoting 
the real squeezing parameter, then the time-evolution of the state is given by 

pit) = U (t)p(0)U j (t). (15) 

Knowing p(t), the Wigner function at time t can be derived by first calculating the 
symmetric characteristic function as the expectation value of the displacement operator 
D(x) as 

e 



Tr 



P(t)D( X ) 



exp [i^/iz cos(u t) — ir]fi x sin(uot)} exp 
2 / cos 2 (wo i sin 2 (wot) 



cos 2 (oj t) ^ sin 2 (c<; t) 



x exp 



' 4 



2*2 



+ 



exp 



sin(2a;ot) 



(16) 

where we substitute e~ 2s = 2<r 2 and e 2s = 2a 2 in the coefficients of £ 2 , rf and ^n. 
Substituting ( IT6|) into the definition of the Wigner function given by ( TTT|) and performing 
two Gaussian integrals over ^ and rj then yields to the SP- Wigner function as 



wS P (x,p;t) 



I exp j - 1 [ r - R(-u t)y\ T ■ T~ l (t) ■ [r - R(-u | 5 



(17) 



where the superscripts "T" and "—1" denote matrix transposition and inversion 
operations, respectively. The suffix SP denotes the Schrodinger picture and the index 
in Wo(x,p;t) represents the fact that the system evolves freely in time. The vectors n 
and r represents the initial displacement of the Wigner function and the coordinates of 
phase space, respectively, as 





;is) 



The matrix T 1 (t) is defined as the clockwise rotation of the squeezing matrix through 
the angle u t 

T-\t) =R(-u t) -E- 1 -R(uot), (19) 
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Figure 1: Free time-evolution of the ^-contour of the SP-Wigner function of an ideal 
squeezed state given by (|T71) in units of uj^t (a x — 1, o~ p — \,n x = 4, in arbitrary units). 
The evolution is a combination of a local counterclockwise rotation of the ellipse about 
its centroid through angle uot, and the global clockwise rotation of its center about the 
origin of the phase space on the dashed-line circle through the same angle. 
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where the diagonal 2x2 squeezing matrix H 1 is given by 

4 

o 4 



(20) 



Notice that the matrix R(u>ot) defines a counterclockwise rotation by uot as 



R(u t)=[ C ° SUJot - s[nUot . (21) 
sin u t cosu t y 

The free-time evolution of the --contour of the SP-Wigner function of an ideal 

e ° 

squeezed state is represented in figure [H As it is illustrated in figures n~atili[ the elliptical 
distribution performs two types of rotations: i) a global clockwise rotation (active 
transformation) of the distribution about the origin of the phase space on a circular 
path defined by the initial amount of displacement; ii) a local counterclockwise rotation 
of the the squeezing direction through the same angle u^t about its centroid. From the 
mathematical point of view, the global rotation is produced by the clockwise rotation 
matrix R(-ojot) acting on the vector /x, whereas the action of the rotation matrix R(uot) 
on the squeezing matrix S" 1 amounts to its local counterclockwise rotations. 

In (fTTj) . if we take R(— oj^t) from left in to the first bracket and R(uot) from right 
in to the second bracket, the resulting Wigner function can be written in terms of a new 
phase space coordinate which is described by 

r = R(u t) r. (22) 

In this new local frame, attached to the centroid of the Wigner function, the local 
counterclockwise rotation of the squeezing direction of the Wigner function about its 
centroid cancels out its global clockwise rotations about the origin 



W HP (xV) = -exp 

7T 



- l - (r' - fxf ■ S- 1 ■ (r - M ) 



(23) 



The passive transformation f l22|) corresponds to changing the reference frame to the one 
in which x(t) and p(t) are rotated back to their initial values x(0) and p(0). The suffix 
HP represents that the Wigner function is described in the Heisenberg picture, in which 
the Wigner function is fixed in time while the operators x and p evolve according to ( 16 oil 
and (l6~fej) . Hence we see that the passive transformation (1221) in phase space corresponds 
to the Heisenberg picture. 

The main message of this section is that the action of the unitary transformation 
Uq on an ideal squeezed state in the Hilbert space corresponds to two types of rotations 
of its Wigner function in the phase space: i) a global clockwise rotation about the origin 
through angle Uot, ii) a local counterclockwise rotation of the Wigner function through 
angle uot about its centroid. Further, the action of Uo on the system in Hilbert space 
is equivalent to a counterclockwise rotation of the coordinate system with angle uot, or 
equivalently, performing a passive transformation in phase space. 
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4. Interaction Picture and mixed rotations 



The interaction picture is the most frequently employed picture in quantum dynamics. 
Specifically, it is the preferred picture in quantum optics where a system is driven by 
the time-dependent electric field of light. 

Considering the single mode harmonic oscillator described in section [2} we now 
assume it undergoes a time-dependent potential V(t). Thus the total Hamiltonian is 
given by 

H = H + V(t) } (24) 

with H being defined in ([3]). V(t) influences the temporal behavior of the system and 
is a linear driving potential described as 

V(t) = g{e~ iullt aa + e iaJl 'a*a t ), (25) 

where g is a real coupling constant. Physically, this corresponds to a forced harmonic 
oscillator with the driving frequency u±. 

In the IP both the state and observables are time-dependent. The observables 
evolve in time according to the Heisenberg equation of motion as 

6 l (t) = ul(t)6u (t), (26) 

where Uo(t) is given by ([TBI . Accordingly, the time- dependent potential in the 
interaction picture, Vr, reads as 

Vj(t) = g( e - int aa + e m a*tf), (27) 

where Q is defined as the sum of the natural frequency uq and the driving frequency oo\: 

n = 000 + 00!. (28) 

When Q = (cox = —oo ), the interaction picture potential becomes time-independent. 
We shall study this case in more details in section 14.21 For the moment we assume 

n ^ o. 

The time evolution of the state is described by the interaction picture unitary 
transformation C/j(£) as 

fr(t) = U^piOplit), (29) 
where Ur(t) is given by the Magnus series [17] 

(30) 



Uj(t) = exp |-i \Axit) + A 2 (t) + A 3 (t) + 
The first three terms in the exponential are given by 



A l (t)= [ dtiV[(ti 
J o 

1 ' dt 2 



2!i 



o 



o 



3!(i 2 ) 



dt 3 / dt 2 
Jo 



(31a) 
(316) 

(31c) 
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where Vr(t) is given by (127)) . Thus the time evolution of the state of the system can 
be derived by calculating Ax(t), A 2 (t), . . ., and subsequently constructing Uj(t), which 
ultimately enables us to calculate p\{t). Then it is straight forward to calculate the 
IP-Wigner function and present a classical correspondence to U\{t). We will ultimately 
transform the IP-density matrix back to the Schrodinger picture pg(t) with the help of 
U(t) as 

p s (t) = U (t) Pl (t)U^t) : (32) 

and study the corresponding SP-Wigner function. 
For the first two term A\(t), and A 2 (t) we have 

ii(t) = i v{t)a ] -v*(t)a + ^ ("*a f - eta) , (33) 

where u(t) is defined as 

u(t) = ~a*e m , (34) 

with a* being the complex conjugate of the displacement parameter defined in (|SJ). Since 
the commutator [a,a^] = 1, the commutator of Vj(t) at two different times turns into a 
purely imaginary function of time. Hence the higher order terms in Magnus series given 
by (131 ap vanish and we are able to write down the exact solution for the equations of 
motion. This conclusion justifies our choice of H\(t) and enables us to illustrate the 
classical correspondence of the IP in phase space. For A 2 {t) we get 

21 12 

A 2 (t) = 9 -^- (SH -sin Qt). (35) 

As a result Uj(t) is given by 

U^t) = e-i^Cnt-asinnt)^ D(2- a *). (36) 

For a system that is initially in an ideal squeezed state described by (TT4|) . the state 
at later time t in the interaction picture reads 

ft(f) = D (y(t)) D(l a *)p(0)D(-l a *)D (-u(t)) , (37) 

where we substituted Eq. (136)) for U\(t) in (1291 . Knowing the state, we can construct 
the IP-Wigner function following the same line of calculations as in section [3] and get 

W}(x,p;t) = -Uxp j-i [r - Ml (t)] T ■ S" 1 • [r - ^(t)] }, (38) 

where the suffix I represents the IP, and the index 1 represents the evolution of the 
Wigner function the takes place in the presence of the interaction Hamiltonian H\(t) 
given by ( 125)) . The vector /Ltj(t) is defined as 
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Figure 2: Time-evolution of the ^-contour of the IP-Wigner function (o~ x = l,a p 

= —2, a = 1,6 = — l,fi = 3,(7 = 5) in one period. The evolution corresponds to 
its global counterclockwise parallel transformation about the origin on the dotted circle 
described by equation ( 1401 . 
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In order to see how the Wigner function evolves in the IP, we look at the equation of 
motion of its centroid given by 

(fo(f)) - /i,) 2 + ((Pl(t)}) 2 = 2 (a 2 + 6 2 )(1 - cosfit). (40) 

As it is shown in figure [2} the time-evolution of the ^-contour of the IP- Wigner function 
corresponds to a displacement of the centroid of its elliptical distribution in phase 
space together with a counterclockwise parallel transformation of the distribution on 
a circle with a time-dependent radius described by equation ( ]4"0~j) . There are two 
important points that can be inferred from figure [2j i) the squeezing direction remains 
invariant during the evolution in the interaction picture, which in our example was 
in the x— direction, ii) The evolution of the IP- Wigner function corresponds to the 
displacement of its centroid around a circle as it was the case for the free Hamiltonian 
evolution of the free evolution of SP- Wigner function of an ideal squeezed state shown 
in figure [TJ The reason is that the IP-unitary transformation Ui(t) is defined as two 
successive displacement operators apart from a phase factor rather than a rotation. It 
is exactly because of this property that for the Hamiltonian (124")) choosing the IP makes 
the study of the time-evolution of the state much simpler. We will see in the next part 
that both of the above mentioned properties get spoiled when we transform back to the 
SP. 



4-1. Back to the Schro dinger picture 

Transforming back to the SP amounts to performing the unitary transformation 

p(t) = U f>i(t)Ul (41) 

which as we showed in section [3] corresponds to a global and a local rotations 
of the corresponding SP-Wigner function through angle uj^t. Then the symmetric 
characteristic function is straightforward to establish. As a result the SP-Wigner 
function of the system in the presence of the interaction Hamiltonian H\(t) is given 
by 



W?* (x, p; t) = I exp | - \ [r - /x s (*)] T ■ T" 1 • [r - /x s (*)] } . 



(42) 



The 2x2 matirx T _1 is given by (fT9|) and the mean value vactor /xg(i) is defined as the 
clockwise rotation of the initial displacement about the origin through angle u^t plus 
two types of transformations of the interaction-induced displacement about the origin: 
a clockwise rotation through angle u^t and a counterclockwise rotation through angle 

UJ\t\ 



/ig(t) = R(-u t) 




R(-u t) , - R(uit) 



.(43) 



The first point we notice by looking at the Wigner function of ( I42p is that when we 
go back to the SP the local and global rotations of the Wigner function appear again 
and this is due to the presence of the unitary transformation U . Furthermore, the 



Visualizing the quantum interaction picture in phase space 



12 




Figure 3: The trajectory of the center of SP-Wigner function in the presence of the 
interaction for different values of u\. ujq is fixed to 1 in arbitrary units. The parameters 
are fixed at a x — 1, a p — |, fi x = —2, a = 1, b = —1, g = 5. 
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global rotation of the centroid of the distribution occurs not on a circle, like in the free 
evolution case illustrated in figure [2j but on a more complex trajectory given by 

((xg(t)> - fi x cosu t) 2 +((p s {t))+fi x sin w t) 2 = 2 (^) 2 (a 2 +6 2 )(l-cos fit). (44) 

Equation (jH]) resembles a glissette [18J. A glissette is defined as the locus of a 
generator point which is moving along a given curve. The generator point in (1441) is 
(fi x cos wot, — l±x sinc^ot) described as the center of the circle in the right hand side of 
( |44l) . This point is rolling along the time-dependent curve given by the left hand side 
of (jHJ). The resulting glissette is shown in figure [3] for different values of Q. 

On top of the rotations induced by Uq, there are two more rotations of the driven 
displacement given by a* due to the presence of the interaction, i) Its clockwise rotation 
through angle u t; ii) its counterclockwise rotation through angle Uit. As a result, the 
SP-Wigner function in the presence of the interaction "dances" in a complicated manner, 
in comparison with its parallel transform in IP, in phase space. This is illustrated in 
figure H] for a fixed value of Q. A comparison between the figure [2] and figure H] makes it 
clear why in studying the evolution of the state of our system of interest the interaction 
picture is preferred. 

Choosing the counterclockwise rotation of the xp-plane given by ( 1221) transforms 
the SP-Wigner function into the IP-Wigner function in the local x'p' frame. This 
corresponds to performing the unitary transformation Uq on the observables in Hilbert 
space. Correspondingly, in phase space replacing r in fT42|) with R(—u t) r' yields the 
same result. Then the vector r — Ms(^) ^ n the local reference frame reads 

r - MS (t) R(-uJot) (r' — fi) — ^(-o; t) ( - b ) + {-b\ (45) 

Substituting f !4^|) into the expression for the SP-Wigner function given by fj4"2j) and 
taking into account that is defined as R(— ujq t)S _1 i?(a;o t), we get 

W^ P (x,p;t) ->• Whx',p';t) = -exp 

7T 



-~ (r' - Ml (t)) T • H- 1 • (r' - Ml (t)) 



(46) 



In fact, by choosing the counterclockwise rotation of the coordinates, we let the 
observables of the system evolve under the free Hamiltonian Hq. Thus the evolution of 
the state in the new local frame is only governed by the interaction Hamiltonian. In 
other words, we separate the time-evolution of the state from that of the observables. 
But this is exactly the definition of the Schrodinger interaction picture |19j . 

We can also transform the Wigner function to the Heisenberg interaction 
picture [19j. Such transformation may be performed by i) a counterclockwise rotation 
of xp frame, i.e. a passive transformation of xj9-frame to a^p'-frame; ii) displacing the 
origin of x'p'-frame by ^( " b ); Hi) counterclockwise rotation of the new origin of the 
x'p'-fra.me with the frequency Q. Thus the transformation from xp-plane to the new 
x"p"-plane is defined as 

r ^ r" = R{u,t)r - ^( * \ + l R (nt) [ \\. (47) 
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Figure 4: Time-evolution of the —contour of W^ p (x,p;t) given by (j42 



4 6 



, [a x = 1, (T p = 

\,Hx = —2, a — l,b — —l,u>i — 2,uj — 1, g — 5). The evolution is a combination 
of local rotations of the center of the ellipse through two different angles uj^t and io\t 
and a global rotation of the ellipse through angle u t. This leads to the "dancing" of 

op 

W? r (x,p;t). The dotted curve describes the motion of the centroid given by 
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Figure 5: A snapshot of the interaction picture Wigner function Wj;(x, p; £)|n=o at some 
time t with the parameters a x = l,a p — h,fi x = —2, a = 1,6= —1,^1 = —&o,g = 5. 
The dotted ellipse represent the Winger function at time t = 0. The dashed line is the 
trajectory of the center of the ellipse. 



The last two transformations are precisely compensating the unitary transformation 
arising from the interaction Hamiltonian. It is then straightforward to describe the 
vector r — /^g (t) in the frame as the clockwise rotations of (r" — fx) through angle 
u t: 

r-^ s {t)^R{-co t){r"-fji). (48) 
Inserting ( T4H1) into ( )42l) gives us the Wigner function in the local x"p" -plane as 



W^ P (x,p;t) ->■ W^ P i 



x ,p ) 



1 



7T 



exp 



.(49) 



We notice that in the local x"p" — frame, only x" and p" are evolving in time according 
to H(t) and the state does not explicitly depend on time. 



4-2. The case Q = 

It is interesting to see what happens when the driving frequency of the interaction 
Hamiltonian given in equation (|25p is equal to that of Hq with opposite sing. The 
first thing we notice is that the interaction picture potential given by equation (|27|) is 
independent of time. As a result the Wigner function is described by 



I 



[x,p;t) 



71 



exp 



1 T — 

- (r-jii(t)|n=o) -S- (r-Aii(*)|n=o) 



(50) 



where the components of /xj(£)|n=o are given by 

(xi(t))\ n=0 = fx x - gbt, (51a) 

(p[(f))|n=o = -</a*. (516) 

Thus the equation of motion for the centroid of the ellipse will be a line instead of a 
circle and the Wigner function will only be parallel displaced along this line in time. 
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The displacement takes place on a line with slope a/b, as it is shown in figure and is 
given by the following equation 

(Pl(«)>|o=o = f «*l(*)>ln=o-M-). (52) 
4-3. Quadratic Hamiltonian 

So far we have considered an interaction Hamiltonian that is linear in the creation and 
annihilation operators. From figures [THS it is apparent that such a Hamiltonian preserves 
the covariance matrix of the Gaussian distribution of the Wigner function, i.e., an initial 
coherent/ squeezed state remains coherent /squeezed under time evolution, although it 
might perform different types of rotations. Mathematically, this means the covariance 
matrix given by (120 j) is invariant under a linear time-dependent interaction term. In 
this section we consider a quadratic interaction term in the Hamiltonian to illustrate 
the deformation of the covariance matrix. We show that in this case the contour of the 
Wigner function initially squeezed in the p-direction expands into a coherent state and 
then becomes squeezed in the x-direction. We choose the quadratic Hamiltonian to be 

H{t) = uotfa + ik [e 2iuJot a 2 - e^*^) 2 ] , (53) 

where k is a real coupling constant. A physical realization of the above Hamiltonian 
may be the interaction of a coherent light beam with a non-linear optical medium. 
In that contains the non-linear susceptibility and the pump field amplitude. 

The group-theoretical approach to the time-evolution of the Wigner function with a 
quadratic potential has been immensely studied [2Q1I2DE2I231I211E51E51EZ!. In what 
follows we use the Lie algebra of SU(1.1)[201 l21~] realized in terms of a and , 

k Q = l -{a)a + \), K + = l -{tf)\ K.= l -a\ (54) 



such that 



K ,K ± 



±K 



2K . (55) 



Thus in terms of these operators we have a Hamiltonian with linear interaction term 

H{t) = 2uj K + ik (e 2ia ">*#_ - e- 2iwot k^j . (56) 

It is straight forward to calculate the unitary time-evolution operator in the interaction 
picture using (1301) . In this case the interaction picture potential Vr(t) becomes 
independent of time and as a result the Magnus series contains only the first term 
Ai(t) given by ( 131 ah . Thus U\it) reads 

U\{t) = e -t(i<—i<+). (57) 

Uj(t) resembles a squeezing operator with a time-dependent parameter, as a 
generalization of the time-dependent displacement operator in (136]) for the linear 
interaction term in terms of a and a) calculated in the previous section. 

Now we can calculate the state of the system at later time t needed for the 
expectation values of the displacement operator. For the sake of simplicity, we assume 
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Figure 6: Time-evolution of the ^-contour of W} (x,p] t) given by (|59l) in a time period 
of 2% (a x = l,o~ p = h,u = 1,k = 0.1). The amount of the squeezing changes 
in time. We see that in the course of time the initially squeezed in the p-direction 
distribution function becomes squeezed in the x-direction due to the presence of the 
quadratic interaction term. The dashed circle represents the vacuum coherent state. 
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Figure 7: Time-evolution of the contour of - value of W^(x,p; t) given by (I6TI) in 
one period of uot (a x — l,a p — h, k = 0.1, Wo = 1). The distribution function starts 
to expand in the x-direction and then squeeze in the p-direction. On top of that it 
also performs a clockwise rotation around the origin. The dashed circle represents the 
vacuum coherent state. 
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the system starts its evolution in the IP from a vacuum squeezed state. This way, 
the centroid of the distribution coincides with the origin of the phase space. Then the 
symmetric characteristic function reads 



Tr 



1 



exp 



e 



[a p e 



+ 



T] 



(5f 



We immediately realize that in the presence of a quadratic interaction the elements 
of the covariance matrix given by fl20|) undergo a scaling (conformal) transformation. 
Furthermore, at a certain time the two elements of the covariance matrix become 
identical thus a squeezed states transforms into a coherent state. Substituting (1551) 
into the definition of the Wigner function given by (1TT1) and performing two Gaussian 
integrals over £ and n gives the IP- Wigner function for the quadratic Hamiltonian as 



w} F (x,p;t) = - exp 

7T 



—r 

2 



(59) 



where the index 2 represents the quadratic interaction term added to the free 
Hamiltonian. 1 (t) is the "breathing" covariance matrix defined as 



g_ 







(60) 



This is illustrated in figure El 

As it is showed in section 14.11 the transformation of the time-dependent density 
matrix back to the SP, amounts to a global clockwise rotation of the contour of the 
SP- Wigner function in phase space. In the case of quadratic interaction term this is 
straight forward to prove and yields to 

W? P (x,p; t) = - exp -\ [r T ■ r-\t) ■ r] , (61) 

where rj x (t) is defined as the clockwise rotated breathing covariance matrix, 

T 2 \t) = R(-u t) ■ E 2 \t) ■ R(u t). (62) 

This is shown in figure H It is worth noticing that since the centroid of the Wigner 
function coincides with the origin of the phase space, the local counterclockwise rotation 
of the elliptical distribution around the centroid is canceled out by the global clockwise 
rotations around the origin as illustrated in figure [71 However, this is not the case for 
an ideal squeezed state as it is illustrated in figure [U 



5. Summary 

We presented a pictorial representation of time evolution of a harmonic oscillator 
driven linearly and quadratically in Schrodinger, Heisenberg, and Interaction picture 
(both Schrodinger interaction picture, in which the evolution of the state is due 
to the interaction potential, and Heisenberg interaction picture, where the the free 
evolution is contained in the state ). This has been done by employing the Wigner- 
Weyl representation to map the density matrices and observables in Hilbert space on 
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distribution functions and variables, correspondingly, in phase space. As an example, we 
took the initial state of the system an ideal squeezed one. We showed that in the presence 
of a linear interaction term the time evolution of the corresponding Wigner function 
in the Schrodinger picture amounts to different types of local and global rotations. 
The interaction picture time-evolution is mapped to a parallel transformation of the 
distribution function on a circle around the origin of the phase space. In transforming 
back to the Schrodinger picture, we showed that the evolution of the Wigner function 
becomes more complicated due to the fact that it performs different types of rotations 
with difference frequencies. However, one can transform from one picture to another by 
performing passive and local transformations of the reference frame. Adding a quadratic 
interaction term to the Hamiltonian amounts to squeezing of the state with a time- 
dependent squeezing parameter. Classical picture of such interaction corresponds to the 
breathing of the Wigner function contours. 
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